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LONG-TERM  GOALS 

The  goal  of  is  to  assess  the  relative  contribution  of  local  dissipation  versus  far-held  radiation  of  internal  waves 
generated  by  the  barotropic  tide  flowing  across  topographically  rough  straits  in  the  Western  Pacific. 

OBJECTIVES 

•  Perform  large  eddy  simulation  (LES)  of  a  bottom  boundary  layer  on  rough  topography  to  quantify  near¬ 
bottom  mixing  as  a  function  of  slope  angle,  stratification  strength  and  other  parameters. 

•  Integrate  the  UCSD  LES  Model  with  the  UNC  Ocean  Model  in  a  joint  effort  with  A.  Scotti  to  cover  a 
realistic  range  of  spatial  and  temporal  processes  deemed  important  to  the  near-held  mixing  observed  in 
the  held  experiment. 

APPROACH 

A  non-hydrostatic  model  that  numerically  solves  the  unsteady,  three-dimensional,  primitive  equations  is  used. 
Advanced  models  such  as  the  dynamic  mixed  model  and  the  dynamic  eddy  viscosity  model  are  utilized  to 
represent  subgrid  processes  in  the  LES  approach.  A  novel  near-wall  model  has  been  developed  so  as  to  increase 
the  Reynolds  number  of  boundary  hows  to  realistically  large  geophysical  values. 

WORK  COMPLETED 

During  the  hrst  year  of  the  grant,  significant  effort  was  expended  in  advancing  the  capabilities  of  the  numerical 
solver.  This  has  allowed  us  to  contrast  the  boundary  layer  among  the  cases  of  subcritical,  critical  and 
supercritical  slope. 

Model  problem 

The  near-bottom  how  resulting  from  a  current  oscillating  with  the  M2  tidal  period  of  12.4  hrs.  on  an  inclined 
surface,  as  illustrated  in  Fig  1,  is  simulated.  The  bottom  is  adiabatic  while  there  is  a  background  thermal 
stratihcation  with  constant  buoyancy  frequency,  N®.  The  how  is  forced 
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Oscillatory  Pressure  Gradient 


by  an  imposed  pressure  gradient  along  the  horizontal  direction  resulting  in  a  barotropic  current 
of  Uosm(Qtd)  in  the  background.  Here  coordinates  x,  y  and  z  denote  the  horizontal  streamwise, 
horizontal  spanwise  and  vertical  directions,  and  u,  v  and  w  denote  the  corresponding  velocity 
components.  The  coordinates  £,  (  and  r)  correspond  to  the  curvilinear  grid  that  conforms  to  the 
bottom  topography. 

The  dimensional  parameters  that  set  the  flow  and  the  thermal  held  are  as  follows:  the 
freestream  velocity  amplitude  f/0,  tidal  frequency  fi,  background  density  gradient  dpbJ  dzd\ooi 
and  the  fluid  properties:  molecular  viscosity  is,  thermal  diffusivity  k,  and  density  p  .  The 
parameters  that  set  the  geometry  are  as  follows:  slope  angle,  (5,  the  slope  length  in  ^-direction, 
and  the  nominal  water  depth,  H .  The  Coriolis  parameter  is  taken  to  be  /  =  0.  There  are 
a  number  of  important  nondimensional  parameters  that  govern  the  how  over  the  slope.  The 
values  of  tidal  excursion  parameter,  Ex  =  Uq/{uI )  and  the  ratio  of  topography  height  to  the 
water  depth  are  taken  to  be  small  values  to  correspond  to  tidal  how  over  a  slope  in  deep  water. 
Moderate  values  of  N  are  considered.  As  a  hrst  step,  the  Reynolds  number  is  small  corresponding 
to  unsteady  transitional  how  rather  than  full-blown  turbulence  at  the  bottom. 

The  angle  9  that  internal  wave  phase  lines  forced  at  the  tidal  frequency,  fl,  makes  with  the 
horizontal  is  given  by 

The  value  of  the  slope  criticality  parameter,  e  =  tan (/3)/ tan(0)  determines  the  slope  response. 
Three  cases  have  been  simulated:  subcritical  e  =  0.5,  near  critical  e  ~  1  and  supercritical 
e  =  1.7. 

Numerical  algorithm  development 

In  the  streamwise  (x)  direction,  the  spatial  discretization  has  been  changed  from  Fourier  col¬ 
location  and  boundary  conditions  have  been  changed  from  periodic,  in  order  to  allow  streamwise 
(x)  development,  key  to  the  boundary  how  on  a  slope.  The  spatial  discretization  corresponds 
to  a  staggered  finite-difference  grid  in  the  x,  z  directions  and  Fourier  collocation  in  the  span- 
wise,  y,  direction.  Spatial  y-derivatives  are  computed  in  spectral  space  while  those  in  the  other 
two  directions  are  computed  using  second-order  finite  difference.  Time  advancement  is  accom¬ 
plished  using  a  third-order  Runge-Kutta  method.  The  viscous  derivatives  are  treated  implicitly 
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using  an  alternating  direction  implicit  (ADI)  method  for  improved  stability  characteristics.  The 
boundary  conditions  on  the  velocity  are  as  follows:  Inflow/outflow  boundary  conditions  along 
with  sponge  regions,  no-slip  at  the  bottom  boundary  for  the  resolved,  moderate-Reynolds  num¬ 
ber  cases  considered  here,  sponge  region  at  the  top  boundary  and  periodic  at  the  lateral,  side 
boundaries.  An  equation  for  the  deviation  from  the  background  density  is  numerically  solved 
with  boundary  conditions  correspond  to  sponge  regions,  periodicity  in  the  lateral  direction  or 
prescribed  slope-normal  gradient  at  the  bottom. 

In  order  to  handle  bottom  topography,  the  Navier-Stokes  equations  have  been  recast  in 
generalized  curvilinear  coordinates.  The  governing  equations  are  transformed  to  a  strong  con¬ 
servation  law.  A  curvilinear  grid-generation  module  has  also  been  added  to  design  a  grid  that 
conforms  to  the  bottom  boundary.  The  grid  is  clustered  in  the  near-bottom  region. 

The  Navier-Stokes  equations  are  solved  using  a  fractional  step  method.  The  pressure  correc¬ 
tion  satisfies  a  Poisson  equation  which,  in  the  present  case,  takes  the  form  of  a  two-dimensional 
Helmholtz  equation.  A  multigrid  solver  is  employed  with  the  following  characteristics:  incom¬ 
plete  L-U  decomposition  as  the  smoother,  weighted  9-point  prolongation  and  restriction,  and 
sawtooth  cycling,  that  is,  smoothing  only  during  grid  transfers  that  occur  during  the  prolonga¬ 
tion  phase  of  the  multigrid  cycle. 

RESULTS 

Near-critical  case,  e  ~  1 

Phase  lines  of  the  internal  waves  at  the  fundamental  frequency,  D,  of  the  M2  tide  are  parallel 
to  the  slope  when  e  ~  1  leading  to  the  phenomenon  of  resonant  generation.  A  tidal  beam 
emanates  to  the  right  as  shown  in  Fig.  2.  Owing  to  the  baroclinic  internal  wave  component, 
the  velocity  in  the  boundary  layer  is  amplified  relative  to  the  peak  barotropic  value  so  that  the 
kinetic  energy  is  as  large  as  ten  times  larger  than  the  background  barotropic  value. 

The  density  held  has  a  response  that  is  asymmetric  between  uphow  and  downflow  as  shown 
in  Figure  3  for  the  near-critical  case.  During  uphow,  there  is  positive  vertical  strain  and  the 
stratification  decreases  while  during  downflow  the  stratification  increases.  The  response  of  the 
disturbance  kinetic  energy  to  such  a  periodic  modulation  of  the  stratification  is  currently  under 
investigation. 

Supercritical  case,  e  =  1.7 

Internal  waves  are  generated  at  the  upper  portion  of  the  slope  where  the  topography  transi¬ 
tions  from  supercritical  slope  to  horizontal.  Similar  to  the  case  with  e  ~  1,  there  is  a  significant 
tidal  beam  associated  with  the  near-critical  extent  of  the  slope  near  the  top  as  shown  in  Fig¬ 
ure  4.  However,  an  important  difference  is  the  presence  of  energetic  beams  at  steeper  angles. 
The  power  spectrum  of  the  streamwise  velocity  in  Figure  5  shows  clear  evidence  of  harmonics, 
nfl  of  the  tidal  frequency.  The  kinetic  energy  beams  with  steeper  angles  correspond  to  these 
harmonics. 

IMPACT  /APPLICATIONS 

Three-dimensional,  highly-resolved  simulations  of  the  bottom  boundary  layer  under  a  barotropic 
tide  on  a  slope  have  been  performed.  The  difference  between  subcritical,  critical  and  supercrit¬ 
ical  slopes  with  respect  to  the  mean  velocity  and  the  internal  wave  held  has  been  highlighted. 
The  impact  of  slope  criticality  on  turbulence  and  dissipation  rates  will  be  obtained  during  the 
second  year. 
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Figure  2:  Resonant  generation  in  the  near-critical  case,  e  ~  1  shown  by  contours  of  E/Ef , 
where  E  is  the  kinetic  energy  whose  freestream  value  is  Ef  =  0.5 Uq.  The  tidal  phase  is  denoted 
by  (j)  and  the  peak  slope-parallel  velocity  is  Up.  The  values  of  </>  for  (a),  (b),  (c)  and  (d)  are 
37t/4,  7 r  (zero  freestream  velocity),  97t/8  and  57t/4,  respectively,  and  corresponding  values  of 
peak  velocity  are  Up/U0  =  2.45,  2.83,  3.62  and  3.01,  respectively. 
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Figure  3:  Top  panel  corresponds  to  a  phase  of  7t/8  where  there  is  upflow  with  peak  velocity 
Up/Uq  =  3.6:  (a)  Kinetic  energy  contours,  (b)  Ispopycnals.  Bottom  panel,  phase  of  97t/8 
corresponds  to  downflow  with  Up/U0  =  —3.62:  (a)  Kinetic  energy  contours,  (b)  Ispopycnals. 
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Figure  4:  Generation  in  a  supercritical  case,  e  =  1.7.  Contours  of  E/Ef,  where  Ef  =  0.5 Uq,  are 
shown.  The  tidal  phase  is  (a)  </>  =  37t/4,  (b)  </>  =  tt,  (c)  0  =  57t/4  and  (d)  </>  =  37t/2,  respectively. 


Figure  5:  Power  spectra  of  the  streamwise  velocity  (semi  log-scale)  as  a  function  of  frequency 
in  the  case  with  e  =  1.7  in  (b).  The  corresponding  spatial  positions  are  shown  by  three  small 
white  horizontal  lines  in  (a). 
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